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This paper reports results of a network theory approach to the study of the United States patent 
system. We model the patent citation network as a discrete time, discrete space stochastic dynamic 
system. From data on more than 2 million patents and their citations, we extract an attractiveness 
function, A(k, I), which determines the likelihood that a patent will be cited. A(k, I) is approximately 
separable into a product of a function Ak(k) and a function Ai(l), where k is the number of citations 
already received (in-degree) and I is the age measured in patent number units. Ai(l) displays a peak 
at low I and a long power law tail, suggesting that some patented technologies have very long-term 
effects. Ak(k) exhibits super-linear preferential attachment. The preferential attachment exponent 
has been increasing since 1991, suggesting that patent citations are increasingly concentrated on 
a relatively small number of patents. The overall average probability that a new patent will be 
cited by a given patent has increased slightly during the same period. We discuss some possible 
implications of our results for patent policy. 

Keywords: innovation, patents, citation network, preferential attachment, aging 



I. INTRODUCTION 



Innovation plays a key role in economic development 
and the patent system is intended (and Constitution- 
ally required in the United States) to promote innovative 
progress. The patent system promotes innovation by giv- 
ing inventors the power to exclude others from using their 
inventions during the patent term. The power to exclude 
is a double-edged sword, however, benefiting the original 
inventor, but imposing costs on later innovators seeking 
to build on past inventions. Thus, the proper design of 
the patent system is an important matter - and a mat- 
ter of considerable current debate. See, e.g., [ElilEl- 
Advances in computer technology and the availability of 
large patent databases have recently made it possible to 
study aspects of the patent system quantitatively. To 
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date the empirical analysis of the patent system has been 
undertaken by economists and some legal scholars. See, 
e -g-j 00- D3 Because patents and the citations between 
them can be conceptualized as a growing network, how- 
ever, techniques from statistical physics that have been 
used in the study of complex networks can be usefully 
applied to the patent citation network 0, Q| . In this pa- 
per we present what we believe to be the first results of 
a network theory approach to the patent system. We ex- 
plore the kinetics of patent citation network growth and 
discuss some possible implications for understanding the 
patent system. 

The paper is organized as follows: In Section II we pro- 
vide background on the United States patent system and 
describe the citation data that is used in this study. In 
Section III we describe a general framework for modeling 
the kinetics of citation networks, define an "attractive- 
ness function" for the evolving network and introduce an 
iterative method for extracting the attractiveness func- 
tion from the data. In Section IV we apply this approach 



2 



to analyze the US patent citation network and explore the 
changes in the kinetics from 1976 to 2000. In Section V 
we discuss some possible implications of our results, and 
mention directions for future research. 



II. PATENTOLOGICAL BACKGROUND 

While a similar approach could be applied to many 
patent systems, including the very important European 
and Japanese patent systems, we begin our analysis with 
the United States patent system for which an extensive 
database of citations has been made available through 
the work of economists Hall, Jaffe, and Trajtenberg pfl. 

An application for a U.S. Patent is filed in the U.S. 
Patent and Trademark Office (USPTO). A patent ex- 
aminer at the USPTO determines whether to grant a 
patent based on a number of criteria, most of impor- 
tant of which for present purposes are the requirements 
of novelty and non-obviousness with respect to existing 
technology Once a patent is issued by the USPTO, it 
is assigned a unique patent identification number. These 
numbers are sequential in order of patent grant. 

Novelty and nonobviousness are evaluated by compar- 
ing the claimed invention to statutorily defined categories 
of "prior art" , consisting in most cases primarily of prior 
patents. Patents are legally effective only for a limited 
term (currently twenty years from the date of applica- 
tion), but remain effective as "prior art" indefinitely. In- 
ventors are required to provide citations to known ref- 
erences that are "material" to patentability, but are not 
required to search for relevant references (though they or 
their patent attorneys often do so) . During consideration 
of the application, patent examiners search for additional 
relevant references. 

Patent citations reflect potential prior art that was 
considered by the examiner. They thus reflect the judg- 
ment of patentees, their attorneys, and the USPTO 
patent examiners as to the prior patents that are most 
closely related to the invention claimed in an application. 
Patent citations thus provide, to some approximation, a 
"map" of the technical relationships between patents in 
the U.S. patent system. This "map" can be represented 
by a directed network, the nodes being the patents and 
the directed edges the citations. Our research uses a 
statistical physics approach inspired by studies of other 
complex networks to attempt to gain insight from that 
"map". 

The patent database we use for the analysis in this pa- 
per was created by Hall, Jaffe and Trajtenbergbased 
on data available from the US Patent Office 0. It 
is available online at http://www.nber.org/patents/ 
The database contains data from over 6 million patents 
granted between July 13, 1836 and December 31, 1999 
but only reflects the citations made by patents after Jan- 
uary 1, 1975: more than 2 million patents and over 16 
million citations. Citations made by earlier patents are 
also available from the Patent Office, but not in an elec- 



tronic format. The Hall, Jaffe and Trajtenberg database 
also contains additional data about the included patents, 
which is described in detail in 0. 



III. MODELING PATENT CITATION 
NETWORKS 

A. Defining the model framework 

In this section we define the mathematical model 
framework we will use for studying patent citations. This 
framework is a discrete time, discrete space stochastic 
dynamic system. Time is measured in patent number 
units. We often "bin" the data from a range of patent 
numbers to obtain sufficient statistics for the analysis. In 
our model, each patent is described by two variables: 

1. k, the number of citations it has received up to the 
current time step and 

2. Z, the age of the patent, which is simply the differ- 
ence between the current time step (as measured 
in patent numbers) and the patent number. Be- 
cause a given patent may cite more than one other 
patent, several citations may be made in one time 
step. 

These two variables define what we call the "attrac- 
tiveness" of a patent, A(k,l) which determines the like- 
lihood that the patent will be cited when the next ci- 
tation is made. In every time step the probability that 
an older patent will be cited is proportional to the older 
patent's attractiveness multiplied by the number of cita- 
tions made in that time step. We find that this simple 
model gives a very good approximation of the observed 
kinetics of the growth of the patent citation network. 

More formally, the state of the system is described by 
ki(t) and h(t), (1 < i < N), where N is the patent 
number of the last patent studied and fcj(i) and l%(t) are 
the in-degree and age, respectively, of patent i at the 
beginning of time step t. The attractiveness of any node 
with in-degree k and age I is denoted by A(k, I). A(k, I) 
is defined such that the probability that node i will be 
cited by a given citation e in time step t is given by 

P[e cites node i] = A ( k ^ l ^\ (1) 

where S(t) is the total attractiveness of the system at 
time step t. 

i 

S(t) = Y,A(k j (t),l j (t)). (2) 
i=i 

The total probability that node i will be cited in time 
step t is thus E(t)A(k i (t),l i (t))/S(t), where E(t) is the 
number of citations made by patent t. A(k, I) and S(t) 



3 



arc defined up to an arbitrary normalization parameter. 
To normalize, we arbitrarily define A(0, 1) = 1. With 
this normalization, S(t) is the inverse probability that a 
"new" node, with k = and 1 = 1, will be cited by a 
given citation during the next time step. 

The A(k, I) function determines the evolution of the 
network, ft describes the average citation preferences 
of the citing patents (the inventors and patent examin- 
ers in reality). In this study, we measure and analyze 
A(k,l) for the United States patent system during the 
time period covered by our data. We find first that the 
parameterization by k and I consistently describes the av- 
erage kinetics of the patent citation network. Of course, 
underlying patent citations are patentee and patent ex- 
aminer evaluations of the significance of the cited patent 
and the technological relationship between the citing and 
cited patents that our probabilistic approach cannot cap- 
ture. The way in which these "microscopic dynamics" 
are translated into the average behavior that we observe 
remains an open question. 

In the following part of this section we will explain our 
method for measuring the A(k, I) and S(t) functions for a 
given network. We believe that this method may be use- 
fully applied to other networks as long as the necessary 
data is available. 



B. Measuring the attractiveness function 

Let us assume that edges are added to the system one 
after another in a fixed order; if two edges are added in 
the same time step (i.e., by the same citing patent), their 
order is fixed arbitrarily for the measurement. Let e be 
an edge and let c e {k,l) be indicator random variables, 
one for each (e,k,l) triple, (1 < e < E tot ,k > 0,1 > 0), 
where E tot is the total number of edges in the system. 
c e (k, I) is one if and only if edge e cites a (fc, I) node (i.e., 
a node having in-degree k and age I) and zero otherwise. 
The probability that edge e cites a (k,l) node, i.e., that 
c e (k,l) is one, is thus given by 



P[c e (k,l) = l] = 



N(t(e),k,l)A(k,l) 

sWf) 



(3) 



where t{e) is the time step during which edge e is added, 
S(t(e)) is the total attractiveness of the system right be- 
fore adding edge e, and N(t(e),k, I) is the number of (k, I) 
nodes in the network right before adding edge e. We thus 
have a formula for A(k,l): 



A(k,l) 



P[c e (k,l) = l]S(t{e)) 
N(t(e),k,l) 



(4) 



In 10} it is easy to determine N((t(e),k,l) for any 
(e, k, I), but S(t(e)) is unknown. Moreover, we have only 
a single experiment for c e (k, I) which is not enough to ap- 
proximate P[c e (k, I) = 1] properly. To proceed further, 
let us define a new set of random variables, each of which 



is a simple transformation of the corresponding c e (k,l) 
variable: 

A -^= c 4mSr - »"«■>.*.<> ><> v> 

If N(t(e), k, I) = then A e (k, I) is not defined. It is easy 
to see that the expected value of any A e (k, I) variable (if 
defined) is A(k, I); thus we can approximate A(k, I) by 



A(k,l) 



1 c e (k,l)S(t(e)) 



E(k,l) ^ N(t(e),k,l) 



(6) 



Here E(k, I) is the number of edges for which 
N((t(e),k,l)) > for any t(e), and c e (k,l) is the real- 
ization of c e (k,l) in the network being studied. 

To calculate this approximation for A(k, I) we need 
to determine S(t(e)), which itself is defined in terms 
of A(k,l). To determine A(k,l) and S(t(e)) self- 
consistcntly, we use the following iterative approach: 

1. First we assume that So(t) is constant, and use © 
to compute Aa{k,l), normalizing the values such 
that A (0, 1) = 1. 

2. Then we calculate Si(t) for each t based on Ao(k, I) 
and use this to determine A\[k, I). 

3. We repeat this procedure until the difference be- 
tween S n (t) and S n +i(t) is smaller than a given 
small e for all t. 

To check this iterative method, we have applied it to 
various well-known models of growing networks, such as 
the Barabasi- Albert model 0] . In these tests the method 
yielded the correct form of the A(k, I) function, which, 
for the BA-model, for example, is A(k, I) = k + a. While 
these tests gave very good agreement overall they also 
suggested that the method cannot accurately measure 
the attractiveness of young nodes (small I) with high in- 
degree (high k), as these occur very rarely in any finite 
sample network. 



IV. RESULTS 

A. The attractiveness function 

The analysis method described in the previous section 
was applied to the patent citation network and the forms 
of S(t) and A(k, I) were determined. Figurcs[I]and|2|sliow 
sections of the A(k,l) function and Figure [3| shows S(t) 
(which is also the inverse of the probability that a "new" 
node with k = and I = I will be cited). For all the 
figures in this paper we have binned the age values into 
300 bins, each containing 7172 patents. Ages and times 
are measured in patent number units. Figures ^ an d [21 
suggest that, for the patent network, the effects of in- 
degrec and age can be separated to good approximation 
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FIG. 1: The measured attractiveness A(k, I) as a function of 
age I for various fixed values of in-degree, k. The bottom 
figure shows only the decreasing tail on log-log scales. 



and that A(k, I) can be written approximately in the form 



A(k,l) = A k (k) ■ Ml). 



(7) 



While this is a reasonable and useful approximation, it is 
also clear that it is only approximately true, e.g., A(0, •) 
decays faster than A(30, •), see the second plot in Figure 
□ 

The measured Ai(l) function for the patent citation 
network has two major features - a peak at approxi- 
mately 200,000 patent numbers and a slowly decaying 
tail. (The very large absolute values of Ai(l) are a result 
of the normalization, A(0, 1) = 1, and are of no indepen- 
dent significance.) The peak at 200,000 patent numbers 
corresponds to a large number of what might be called 
"ordinary", relatively short-term citations. In 1998- 
1999, 200,000 patent numbers corresponded to about 15 
months. The tail is best described by a power-law decay: 
At (I) ~ r' 3 with (3 w 1.6. The observation of this power 
law decay is an important result. It indicates that while 
typical citations are relatively short-term, there are a sig- 
nificant number of citations that occur after very long 
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FIG. 2: The measured attractiveness A(k, I) as a function of 
in-degree, k, for various fixed values of age, I. 



delays. Very old patents are cited, suggesting that the 
temporal reach of some innovations, which perhaps can 
be described roughly as "pioneer" , is very long indeed. 
Moreover, because Ai(l) is approximately independent 
of k - i.e., approximately the same power law decay is 
observed even for small k - the power law tail of Ai(l) 
demonstrates that there is a significant possibility that 
patents that have gone virtually un-cited for long periods 
of time will reemerge to garner citations. This slow power 
law decay of Ai(l) thus suggests the unpredictability of 
innovative progress. 

The measured A k (k) function increases monotonically 
with k, as Figure |21 suggests. Higher in-degree always 
means higher attractiveness. Since the citation proba- 
bility is proportional to the attractiveness, this means 
that the well-known preferential attachment, or "rich get 
richer" effect is at work here - the more citations a patent 
has received, the more likely it is to receive another. The 
functional form of Ak(k) is a power law over the entire 
range of k values. Ak(k) ~ k a + a, where a = 1.19 
and a = 1.11. We estimated these parameters using the 
smaller values of k, for which we have more data, and 
the first 100 age bins. We then checked the results by 
comparing with more extensive fits. 

Preferential attachment and its variations are well 
studied, see the reviews by Albert and Barabasi [| and 
by Newman [l4|. Linear preferential attachment (a = 1) 
without aging has been shown to result in a degree distri- 
bution (frequency of nodes with degree k) with a power 
law tail Krapivsky et al. H3 have studied nonlin- 

ear preferential attachment. In the model they studied 
there was no aging, A(k, I) = A k (k) = k a + a. For a > 1, 
as is observed in the patent citation network, their cal- 
culations predict a condensation of node connectivity, in 
the sense that with high probability most of the edges are 
connected to only a small number of nodes. More specif- 
ically, in their model, if (m + \)/m < a < m/(m — 1) the 
number of nodes with more than m incoming edges is fi- 
nite, even in an infinite network. For the patent network 



5 



1976 1980 1984 1988 1991 1994 1997 1999 




dependence: 



4000000 



4500000 5000000 5500000 6000000 
time [patent number] 



FIG. 3: The total attractiveness S(t) of the patent network 
versus time in units of patent numbers. For ease of reference 
the time in years is indicated by filled circles and vertical 
lines. The left and right figures in the inset show the total 
attractiveness of Model 1 and Model 2, as discussed in the text 
surrounding equation (8) and (9). The left plot shows results 
of simulations for A = 1/50 (circles) and 1/5 (triangles). The 
right plot shows results of simulations for n — 1/2 (circles) 
and 1 (triangles). The simulations were run for 10,000 nodes. 



7/6<a<6/5 suggesting that, if there were no aging, the 
number of patents receiving more than 6 citations would 
be very small, though those patents would account for a 
large fraction of all of the citations. Aging complicates 
this picture, of course, and likely precludes a complete 
condensation onto a few nodes. However, the fact that 
the observed preferential attachment is super-linear does 
indicate a tendency toward what might loosely be called 
"stratification" - many nodes with very few citations and 
a few nodes with many citations. 



B. The total attractiveness 

The total attractiveness function, S(t), (see Fig. [3J| of 
the US patent system increases with time. The initial 
steep increase is only a finite size effect and comes from 
the fact that the citations made by pre-1975 patents are 
missing from our database. From about 1984 on, how- 
ever, S(t) displays a slow but steady increase. One way 
to interpret this increase is that the probability that a 
patent will be cited by a given citation (which is propor- 
tional to 1/S(t)) is decreasing as the size of the network 
increases. The decrease is determined in part by the rate 
at which patents age, which determines the number of 
patents "available" for citation. 

To better understand the behavior of the S(t) func- 
tion, we simulated two simple growing network models 
with two different "toy" A(k, I) functions, with linear 
preferential attachment and two different forms of age 



Model 1 
Model 2 



A^kJ) = (fc + l)-e" A ' and 
A 2 (k,l) = (fc + !)•/-". 



(8) 
(9) 



Si(t) and 5*2 (t) were determined for these models, see 
Fig. |3 inset. In these models a single edge was added to 
the network at each time step. 

When the attractiveness function decays exponentially 
with age (Model 1), the total attractiveness fluctuates 
around a constant value, which is determined by A and is 
independent of the system size. In Model 1 the probabil- 
ity that a new node will be cited by a particular citation 
is thus always the same. The exponentially decaying age 
dependence means that the effect of very old nodes is 
negligible; there is effectively a constant-sized "band" of 
recent nodes that remain "citable" . In a patent citation 
system, such an exponential age dependence would im- 
ply, contrary to our observations, that the importance of 
innovations is short-lived. 

In Model 2, the behavior of S2OO depends on [i. If 
\x is below a limit value [i\ (which is about 1), Si(t) is 
sharply increasing. If \i is between \i\ and another limit, 
fi2 (which is about 1.5) S^C*) increases slowly. When 
li is higher than /Z2, ^(t) fluctuates around a constant 
value, as does S\(t). Thus Model 2 exhibits a crossover 
from a regime of slowly decaying age dependence in which 
old nodes remain influential, to a regime of more rapidly 
decaying age dependence in which old nodes are "forgot- 
ten" . These results are in good agreement with other 
theoretical studies about aging and preferential attach- 
ment; see the work by Dorogovtsev and Mendes £|, Zhu 
et al. 0, Klemm and Eguluz [Toj j. 

These simple models would suggest that the patent ci- 
tation network (which has an aging exponent of about 
1.6 - above /Z2 - should have S{t) roughly constant in 
time. Indeed, a third toy model with superlinear pref- 
erential attachment (exponent 1.2) and power-law aging 
(exponent 1.6) displayed a roughly constant S(t). How- 
ever, the observed S(t) for the patent system increases 
with time. 

The explanation for this apparent contradiction lies 
in the fact that the average number of citations made by 
each patent (and hence, since we measure time in units of 
patents, the number of citations made in each time step, 
E(t)), has increased approximately linearly with time in 
the real patent citation network, e.g., it was 4.69 in 1975 
and 10.66 in 1999. See Fig. 0] A modified simulation 
with superlinear preferential attachment (exponent 1.2) 
and power-law aging (exponent 1.6), but with linearly 
increasing number of citations per time step confirmed 
that, as in the real patent data, the increasing number of 
citations made in each time step results in an increasing 
overall attractiveness, S(t). 

The probability that patent i will be cited in a given 
time step (in other words, by a particular patent rather 
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FIG. 4: The top figure shows the number of citations made 
per patent, E{i), as a function of time in units of patent 
number. The bottom figure shows E(t) / S(t), corresponding 
also to the probability that a new patent with k — 0, 1 = 1 
will be cited, as a function of time in units of patent number. 



than by a particular citation) is 



P[ki(t+l) = ki{t) + l]=E(t) 



A(h(t),h(t)) 
S(t) 



(10) 



The probability that a new patent (k = 0, 1 = 1) will be 
cited by the next patent is thus given by E(t)/S(t), which 
is shown in Fig. 21 From this plot one can see that the 
increase in the number of citations being made outweighs 
the increase in S(t), so that the probability that a new 
patent will be cited has increased over time, despite the 
increasing S(t). Patents do not get "lost in the crowd" 
the way wc might have predicted from the simple models. 
Instead, patentees and patent examiners have on average 
increased the number of citations made by each patent 
to more than compensate for the increasing S(t). 



C. Change in the patent system dynamics 

While it is well known that there has been a signifi- 
cant increase in the number of US patents granted per 
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FIG. 5: The measured value of a as a function of time, mea- 
sured as described in the text. The time in years is indicated 
as in Figure 2] 



year since 1984 |6|,|9j, the underlying reason for this in- 
crease is not clear. Has there simply been an acceleration 
of technological development in the last twenty years or 
has there been a more fundamental change in the patent 
system, perhaps, as many have suggested, as a result of 
increased leniency in the legal standard for obtaining a 
patent [i| . A complete answer to this question is far be- 
yond the scope of the present investigation. However, 
our kinetic model does permit us to ask whether there 
has been any deep change in the growth kinetics of the 
patent citation network. Because we measure time in 
units of patent number, a mere acceleration of techno- 
logical progress should leave A(k, I) unchanged in patent 
number "time". A change in A(k, I) indicates some other 
source of change. 

Thus far, we have assumed a time-independent A(k, I), 
which is reasonably consistent with our observations. In 
this section, we relax this assumption to ask the more 
subtle question of whether there has been a change in 
patent system kinetics over and above the acceleration 
that is already reflected in our choice of time units. 
Specifically, wc allow a and (3 to vary with time and 
ask whether there has been a significant change in these 
parameters between 1980 and 2000. 

To answer this question we measured the parameters 
of the system as functions of time. To perform the 
fits, we averaged over a 500,000-patcnt sliding time win- 
dow and calculated the parameters after every 100,000 
patents. The measured a parameters arc plotted in Fig- 
ure[5| There is a significant variation over time. The time 
dependence of the important (3 parameter was also ex- 
plored, but no significant time dependence was observed 
to within the statistical errors. 

The plot for the a parameter shows that there are two 
regimes. In the first regime, prior to about 1991, a is de- 
creasing with time, while in the second, starting around 
1991, there is a significant increase. As noted earlier, 
the a parameter has some very important consequences 
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for the growth of the network: the higher a, the more 
"condensed" or "stratified" the network will be. The 
increasing a in the patent citation network indicates in- 
creasing stratification - a smaller and smaller fraction of 
the patents is receiving a larger and larger fraction of the 
citations. This change is not simply a result of acceler- 
ating numbers of patents being granted, but suggests a 
more fundamental change in the distribution of patents 
that are being issued. 

V. CONCLUSIONS 

We have presented a stochastic kinetic model for 
patent citation networks. Though a complex process un- 
derlies each decision by a patent applicant or examiner 
to cite a particular patent, the average citation behavior 
takes a surprisingly simple form. The citation probabil- 
ity can be approximated quite well by the ratio of an 
"attractiveness function", A(k,l), which depends on the 
in-degree, k, and age in patent numbers, I, of the cited 
patent, and a time-dependent normalization factor, S(t), 
which is independent of k and I. 

We introduced a method to extract the A(k, I) and S(t) 
functions of a growing network from a specification of 
the connection history. We applied this technique to the 
patent citation network and, though no assumptions were 
made as to the functional form of A(k, I), the measured 
A(k, I) function was well described by two approximately 
separable processes: preferential attachment as a func- 
tion of in-degree, k, and power law age dependence. The 
interplay of these two processes, along with the growth 
in the number of citations made by each patent, governs 
the emerging structure of the network. Particularly note- 
worthy are our finding that the preferential attachment 
is super-linear, implying that patents are highly strati- 
fied in "citability" , and our finding of a power law tail in 
the age dependence even for small k, indicating not only 
that some patents remain important for very long times, 
but also that even "dormant" patents can re-emerge as 
important after long delays. 

We also used our technique to investigate the time de- 
pendence of the growth kinetics of the patent citation 
network. Overall, we find that the increasing number 
of patents issued has been matched by increasing cita- 
tions made by each patent, so that the chance that a 



new patent will be cited in the next time period has even 
increased over time. This result suggests that on average 
patents arc not becoming less "citable". However, we 
also find that there has been a change in the underlying 
growth kinetics since 1991. Since 1991, preferential at- 
tachment in the patent system has become increasingly 
strong, indicating that patents are more and more strat- 
ified, with fewer and fewer of the patents receiving more 
and more of the citations. A few very important, per- 
haps "pioneer", patents seem to dominate the citations. 
This trend may be consistent with fears of an increas- 
ing patent "thicket" , in which more and more patents 
are issued on minor technical advances in any given area. 
These technically dense patents must be cited by patents 
that build upon or distinguish them directly, thus requir- 
ing that more citations be made, but few of them will be 
of sufficient significance to merit citation by any but the 
most closely related patents. Further work will be needed 
to understand this change in citation network kinetics. 

This work is only the beginning. There are many fur- 
ther applications of network analysis to the patent cita- 
tion network that are likely to bear fruit. It will be pos- 
sible, for example, to compare the structural and kinetic 
behavior of the network for patents in different techno- 
logical areas, to investigate the degree of relatedness be- 
tween patents in seemingly disparate technologies, and to 
explore more detailed structural indicators, such as clus- 
tering coefficients and correlation functions. Also, it may 
be possible to compare the growth of patent systems in- 
ternationally perhaps providing a means to distinguish 
between the effects of global technological change and 
those of nation-specific legal changes. Finally, it will be 
interesting to compare the behavior of the patent citation 
network with that of other networks (such as the scien- 
tific journal citations discussed in [lf|) to gain deeper 
insight into the behavior of complex networks in general. 
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